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Abstract. We include the effect of evaporation in our evolutionary calculations of close-in giant planets, based 
on a recent model for thermal evaporation taking into account the XUV flux of the parent star (Lammer et al. 
2003). Our analysis leads to the existence of a critical mass for a given orbital distance m cr i t (a) below which the 
evaporation timescale becomes shorter than the thermal timescale of the planet. For planets with initial masses 
below m cr i t , evaporation leads to a rapid expansion of the outer layers and of the total planetary radius, speeding 
up the evaporation process. Consequently, the planet does not survive as long as estimated by a simple application 
of mass loss rates without following consistently its evolution. We find out that the transit planet HD 209458b 
might be in such a dramatic phase, although with an extremely small probability. As a consequence, we predict 
that, after a certain time, only planets above a value m cr it(a) should be present at an orbital distance a of a star. 
For planets with initial masses above m c tit, evaporation does not affect the evolution of the radius with time. 
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^ ■ 1. Introduction temperature of the planet. The first attempt of L03 to 

model such a complex process yields an escape rate in 
The increasing number of discovered giant planets orbit- good a g re ement with the observational determinations of 
C3 ; ing at < 0.1 AU from their parent star raises fundamental yidal-Madjar et al. (2003) for HD209458b, providing en- 
questions about their formation and migration process and C0U raging support for further exploration. Moreover, since 
about the influence of the parent star through irradiation stellar xuv nuxes vary significantly with time and can 
or tidal effects. The recent discovery of an extended at- be order of magnitudes larger at very young ages, these 
mosphere for the transiting exoplanet HD209458b (Vidal- evaporation rates are much larger at the planet early evo- 
Madjar et al. 2003) highlights the occurence of atmo- mt ionary stages. L03 thus suggest that mass loss could be 
spheric evaporation for these close-in planets. Whether an important event in the life of close-in exoplanets, con- 
such evaporation due to heating from the incident stellar trarily to what was pre viously thought. It is the purpose 
flux leads to major mass loss during the planet lifetime, of this letter to explore this issue by taking into account 
and whether this process affects significantly the struc- consistently the thermal escape rates of L03 along the evo- 
ture of the planet and thus its m-R relationship is an lution of strongly i rra di a ted planets (Baraffe et al. 2003, 
open question, which is of prime importance for our under- hereafter B03) . Since an important issue of this analysis is 
standing of planetary system formation. New evaluations to determine whether evaporation affects significantly the 
of atmospheric thermal evaporation rates by Lammer et inner structure an d thus the m-R relationship of exoplan- 
al. (2003, L03), based on exospheric heating by stellar etS) we focus on the case of preS ently detected transits, 
XUV radiation, yield significantly larger rates than the namely HD209458b, with a = 0.046AU (Charbonneau et 
previous estimates assuming Jeans escape at the effective al 2 000) and OGLE-TR-56b, with a = 0.023AU (Konacki 
et al. 2003). 
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2. Description of the models 

The evolutionary calculations are based on the consis- 
tent coupling between the irradiated atmospheric and in- 
terior structures as described in B03, and in Barman et al. 
(2001) for the atmosphere model calculations. Such a con- 
sistent treatment of the irradiated atmospheric structure 
and the internal, partially radiative structure successfully 
reproduces the observed parameters of the transit planet 
OGLE-TR-56b (Chabrier et al. 2004). 

Details of the model used to derive thermal evapora- 
tion rates can be found in L03. The basic idea relies on 
the fact that the energy deposition by stellar XUV leads to 
exospheric temperatures higher than the blow-off temper- 
ature for H. Therefore, the classical Jeans' description of 
thermal escape must be replaced by a hydrodynamic mod- 
eling of the expansion and mass loss. The energy-limited 
atmospheric mass loss rate M can be written: 



M = 3/3 3 ^/(Gp), 



(1) 



where /? is the ratio between the expansion radius Ri, 
where the bulk of the XUV radiation is absorbed, and the 
planetary radius i?p and p = (3rn)/(47ri?p) is the mean 
planet density. L03 estimate /3 by applying the hydrody- 
namic model of Watson et al. (1981). The term is the 
stellar flux, averaged over the whole planet surface, tak- 
ing into account both the contribution in the 1-1000 A 
wavelength interval and the Lyman-a flux, so that the to- 
tal contribution F+ for an orbital separation a (in AU) is 
F* = (Fxuv + F a )/a 2 . Although (3 is expected to vary 
with F* and a, we fix it at the maximum value (3 — 3 
found in L03, which is in good agreement with the ob- 
served expanded exosphere of HD209458b (Vidal-Madjar 
et al. 2003). Estimates of XUV and Lyman-a fluxes and 
their time dependence are based on the current solar value 
and a collection of data for solar type stars covering an age 
from 100 Myr to 8 Gyr. As in L03, we adopt for the XUV 
and Lyman-a contributions: 



Fxuv(£) 
Fxuv (t) 



6.13^ 119 /xuv 
Fxuv(O.l) 



if t ^ 0.1 Gyr, 
if t < 0.1 Gyr 



F a (t) =3.17i- 075 / Q 
F a (t) =F a (0.1) 



if t ^ 0.1 Gyr, 
if t < 0.1 Gyr 



(2) 



(3) 



where t is the age in Gyr, and /xuv = 8.5 10~ 4 W.m -2 
and f a = 1.42 10~ 3 W.m~ 2 are respectively the Sun XUV 
and Lyman-a surface-averaged fluxes at 1AU (Woods & 
Rottman 2002). Equations (2) and (3) recover the solar 
values for £=4.5 Gyr. 

3. Evolution of evaporating planets 

We have implemented the evaporation rates given by 
Equ. (1) in our evolutionary code for irradiated planets, 
focusing on the two aforementioned orbital separations, 



a=0.023 AU and a = 0.046 AU. We investigate plan- 
ets with initial masses 0.5 Mj to 5 Mj, orbiting a solar 
type star (T e // 4 = 6000 K). For this mass range, evapo- 
ration rates vary from ~ 10~ 8 Mj/yr at young ages to ~ 
10- 12 Mj/yr for t > 5 Gyr. 

3.1. Response to mass loss: the case of a 1 Mj planet 

Figure \\\ shows the evolution of an evaporating planet 
with initial mass lMj. A key quantity to understand 
the response of an object to mass loss is the ratio of 
the mass loss rate timescale tvj = m/M to the thermal 
timescale, characterised by the Kelvin-Helmotz timescale 
£kh ~ 2Gm 2 / (RL). Such a ratio is displayed in the upper 
panel of Fig. QJ As long as t^/tun > 1, the evolution of 
the planet radius R(t) is barely affected by evaporation, as 
seen by the comparison between the non-evaporating and 
evaporating curves in the lower panel of Fig. \T\ (see also 
Fig- El- However, when becomes shorter than £kh, the 
evolution of the planet changes drastically: its mass de- 
creases faster than its contracting radius. Consequently, 
m/Rp decreases increasingly rapidly with time, or con- 
versely the planet mean density increases more and more 
slowly, compared to a case with slower evaporation. The 
evaporation rate thus decreases more slowly with time 
since M oc 1/p oc Rp/m (Eq. (1)). Consequently, the ratio 
%l / £kh keeps decreasing. When t ^ < £kh / 10 , the response 
of the planet becomes rather violent, as indicated by the 
sudden increase of the radius at t ~ 20 Myr for the case a 
= 0.026 AU (Fig. [3 long-dashed curve). This stems from 
the fact that the outer layers start to expand significantly 
and the radius, instead of decreasing slowly or remaining 
constant, increases rapidly, yielding a sudden increase of 
the evaporation rate oc R Pl producing in turn further ex- 
pansion of the outer layers. The situation turns into kind 
of a runaway behaviour and points to a catastrophic fate 
for the planet, at least concerning its hydrogen-rich enve- 
lope. 

The expansion of the outer layers, yielding eventually 
such a violent reaction, can be understood in terms of 
entropy balance, in analogy with mass loss of low mass 
secondaries in compact binary systems (Ritter 1996). The 
entropy profile of an irradiated planet of mass m is con- 
stant throughout most of the structure. The outer lay- 
ers, however, are radiative, and are characterized by a 
nearly isothermal, high entropy profile (Guillot et al. 1996; 
Barman et al. 2001; B03; Chabrier et al. 2004). The mass 
enclosed in the radiative layers is typically ~ 10~ 5 Mj. 
With an evaporation rate M ~ 10~ 8 Mj/yr, it takes 
At ~ 10 3 yr to evaporate all these layers. The upper con- 
vective zones of entropy S conv are then exposed to irra- 
diation of the parent star, and become radiative with a 
significantly higher entropy SVad- In terms of gravother- 
mal energy rate defined as: 



Egrav - Af 



-T 



( S rad $ conv ) 

At : 



(4) 



Baraffe et al.: The effect of evaporation on the evolution of close-in giant planets 



3 



fgrav is thus negative. Gravothermal energy is thus re- 
leased in the outer layers, and converted mainly into ex- 
pansion work. At some time, this expansion work be- 
comes significant compared with the (positive) gravother- 
mal energy rate due to the planet contraction and cooling. 
This occurs when £kh S> tv,, yielding a runaway situation 
where evaporation is speeding up and leading eventually 
to a catastrophic event if other processes do not regulate 
the increase of R and M. When a planet reaches such a 
critical regime, evaporation is speeding up drastically, at a 
much higher rate than expected from standard scenarios. 

As seen in Figure < £kh after 10 Myr for a = 

0.023 AU and after 50 Myr for a = 0.046, respectively. The 
a = 0.046 sequence shown in Fig. Estops before showing 
a strong increase of R simply because it reaches the lower 
Teff-limit of our irradiated atmosphere grid, i.e. T e s = 50 
K. The planet at a = 0.046 AU survives ~ 1 Gyr and 
the one at a = 0.023 AU only 20 Myr. Using the same 
evaporation rates without considering the effects on the 
evolution yields a time for complete evaporation of a 1 
Mj planet of - 2 Gyr at a = 0.046 and ~ 40 Myr at a 
= 0.023 AU. A consistent treatment of evaporation and 
evolution, yielding the aforedescribed evaporation speed- 
ing up process, thus decreases appreciably, by a factor of 
~ 2, the timescale of a planet for complete evaporation. 

3.2. A critical mass for speeding up evaporation 

Following the results of the previous section, we define a 
critical mass m cr it below which — ^kh/10 is reached 
in <5 Gyr. A planet with initial mass m < TO c rit orbit- 
ing a Sun-type star will thus evaporate entirely (or leave 
at least a rocky core) in less than 5 Gyr. Based on the 
present evaporation rates, we find that m cr it ~ 1.5Mj 
for a = 0.046 AU and m crit ~ 2.7Mj for a = 0.023 AU, 
as illustrated in Figure |2 Note that in our calculations, 
evaporation starts at the very beginning of our evolution- 
ary sequences. This is questionable, since heating from the 
parent star is likely to be smaller at the very early stages 
of the planet evolution, because either of the presence of a 
protoplanetary disk or the fact that the planet may form 
at a larger orbital distance and migrate inwards. Both 
processes have timescales obviously linked to the disk life- 
time, with a reasonable upper limit of ~ 10 Myr (Strom et 
al. 1993, Armitage et al. 2003). On such short timescales, 
evaporation effects on the evolution of the planet are un- 
significant in the presently considered mass range. Indeed, 
we checked that starting evaporation at an age f ~ 10 Myr 
does not change the fate of the planet and thus rn C rit- A 
prediction of the present calculations is thus the absence of 
exoplanets below the aforementioned critical masses and 
orbital distances after 5 Gyr. 

3.3. Evolution of the radius 

For planets with initial mass m > m cr j t , the evolution 
of the radius, R(t), is barely distinguishable from non- 




6 7 8 9 10 

log t 

Fig. 1. Evolution of an evaporating planet with initial 
mass 1 Mj at different orbital separations. From lower 
to upper panel: the planet radius (in units of Jupiter ra- 
dius Rj = 7.1492 10 9 cm), mass, evaporation rate M in 
Mj /yr and the ratio of the evaporation rate to thermal 
timescale log (t^/tiai)- The solid (with evaporation) and 
dotted (no evaporation, lower panels) curves correspond 
to a = 0.046 AU. The long-dashed (with evaporation) and 
dash-dotted (no evaporation, lower panels) curves corre- 
spond to a = 0.023 AU. 

evaporating cases. This is illustrated in Figure |21 which 
portrays the evolution of a 2.7 Mj planet located at o = 
0.023 AU from its parent star. With the present evapora- 
tion rates, the evaporating sequence (solid line) reaches 
a mass 1.5 Mj and a radius 1.12 i?j after 3 Gyr, in 
agreement with the observed properties of OGLE-TR-56B 
(Torres et al. 2003). The dashed curve shows the evolu- 
tion of a non-evaporating planet of mass 1.5 Mj, which 
reaches a similar radius 1.11 Rj at t = 3 Gyr. At such 
an age, the evaporation rate is ~ 5.5 10 12 g s _1 , orders 
of magnitude larger than recent estimates of ~ 10 3 gs" 1 , 
based on Jeans escape rates at exobase temperatures close 
to T efi (Sasselov 2003). 

Concerning the case of HD209458b, we cannot find an 
evaporating sequence reproducing its properties. As for 
non-evaporating models (see B03), the predicted radius 
for the observed mass and inferred age of the system is 
about 25 % smaller than the observed value. Interestingly 
enough, we find that a planet at a =0.046 AU with initial 
mass ~ 1.1-1.2 Mj reaches the critical regime i^/iicH < 
0.1 precisely at the age of HD209458b (log t = 9.6 - 9.85) 
(see Fig. 0). Since our analysis suggests a violent reac- 
tion of a planet when reaching this regime, with rapid 
expansion of the outer layers (see Fig. 1), we may won- 
der whether HD209458b is not in such a regime. Although 
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Fig. 2. Ratio of the mass loss timescale to thermal 
timescale t ^ /£kh as a function of time (in yr) for planets 
with different initial masses, as indicated on the panels. 
Upper panel corresponds to orbital separation a — 0.046 
AU and lower panel to a = 0.023 AU. In both panels, 
the solid curve indicates the critical mass below which 
^mAkh < 1/10 in less than 5 Gyr (see text). 

the probability to see a transit planet in this rapid phase 
is very small, the discovery of other similar systems with 
unexplained large radii would suggest further attention to 
this scenario. 



4. Conclusion 

The present calculations, including a consistent treatment 
of irradiation and evaporation due to XUV/Lyman-a irra- 
diation during the evolution of irradiated planets, suggest 
the existence of a critical mass m C ritj which varies with or- 
bital separations, below which the evaporation timescale 
becomes significantly shorter than the thermal timescale 
of the planet. Based on the present evaporation rates for 
solar conditions, we get m C rit ~ 1-5 Mj at a = 0.046 AU 
and m C rit ~ 2.7Mj at a = 0.023 AU. For objects with mass 
below this critical mass, we find that: (i) the planet will 
evaporate entirely (except possibly for the central rocky 
core) within about 5 Gyr; (ii) its lifetime is shortened by a 
factor ~ 2 compared with the time predicted for complete 
evaporation but omitting the effect on the evolution; (iii) 
its outer layers expand rapidly and its radius eventually 
increases at some time, (iv) the planet undergoes a phase 
of rapid mass loss which could expel part or all of its re- 
maining hydrogen-rich envelope in a very short timescale. 
Planets with initial masses m\ > m cr i t survive to evapora- 
tion on a lifetime longer than 5 Gyr, and evolution is simi- 
lar to the case of a non-evaporating irradiated planet. The 




Fig. 3. Effect of evaporation on the radius and mass as 
a function of time (in yr) for planets at a — 0.023 AU 
from their parent star. The solid curve is an evaporating 
planet with initial mass 2.7 Mj, which reaches 1.5 -Mj in 3 
Gyr, reproducing the properties of OGLE-TR-56b (Torres 
et al. 2003). Comparison is made with non evaporating 
evolutionary sequences with mass 2.7 Mj ( dash-dot) and 
1.5 Mj (dash). Note that the solid and dash-dot curves 
are undistinguishable in the lower panel. 

values of m cr ;t depend on the evaporation rates which are 
still very uncertain and rely on a rather primitive model 
for such a complex process. Moreover, applying evapora- 
tion rates determined at the exosphere to the photospheric 
surface of the planet, the relevant boundary for evolution, 
is an extreme simplification. Uncertainties on the mass loss 
rates, however, do not affect the qualitative existence of 
a planet critical mass, for a given orbital distance, below 
which a planet eventually will react violently to evapora- 
tion and will expand again, possibly until complete evap- 
oration. Such a behavior bears important consequences 
on the lifetime of close-in planets and on our understand- 
ing of their mass-orbital period distribution. Our results 
thus provide an excellent motivation to pursue efforts to 
understand evaporation effects and to explore further the 
influence of the high energy part of the parent star inci- 
dent flux. 
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